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diction from  tests  appmrs  to  be  that  the  stress  at  the  half-ligament  di- 

mension, calculated  from  J(^  and  the  Hutchinson-Rlce-Rosengren  singularity, 
be  no  more  than  1.05  to  1.1  times  the  yield  strength.  For  smaller  specimens, 
the  J-lntegral  varied  linearly  with  ligament  size  and  was  surprisingly  similar 
for  single  and  double  grooved  specimens.  On  the  other  hand,  fully  plastic 
trlaxiallty  was  shown  to  be  very  Important  in  welded  plates  of  A-36  steel,  — ^ 
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In  which  the  transition  temperature  for  brittle  crack  propagation  was 
raised  40  K above  the  Charpy  value  by  the  triaxiality  due  to  simulated 
cracks  on  both  sides  of  the  plate.  This  may  also  be  of  interest  in 
f ragmentation 

When  the  J- integral  and  HRR  fields  were  no  longer  valid,  more 
crack  tip  opening  in  terms  of  dimple  diameter,  more  sliding-off,  and 
less  delamination  occurred  in  the  single  grooved  specimens  than  in  the 
double  grooved  ones.  The  dominant  mechanism  of  crack  growth  seems  to 
be  interactions  of  holes  with  the  crack  tip  one  at  a time. 

The  crack  toughness  of  a variety  of  steel  and  aluminum  alloys 
seemed  to  be  more  closely  correlated  with  the  yield  strain  Y/E  than 
with  any  other  micros true tural  or  metallurgical  features  except  cleavage 
fracture.  This  fact,  along  with  a review  of  fracture  mechanisms,  in- 
dicates that  the  primary  possibilities  for  improvement  in  fracture  re- 
sistance are  a)  for  Initiation,  increasing  the  adherence  between  the 
matrix  and  second  phase  particles,  and  b)  for  propagation,  second  phase 
particles  oriented  to  promote  crack  meandering  and  branching. 

The  difficulty  with  analysis  in  the  elastic-plastic  to  fully 
plastic  regimes  suggests  that  more  emphasis  should  be  put  on  realistic 
tests,  such  as  ones  involving  plane  strain  ductility  and  cracked  plates 
simulating  various  weld  defects. 

The  boundary  Integral  method  using  dislocations  in  internal  ele- 
ments for  numerical  calculations  of  plasticity  appears  to  be  useful 
primarily  for  small  plastic  zones  around  cracks  in  large  bodies.  Insta- 
bility of  the  method  eventually  developed  with  non-hardening  plasticity. 
For  rate-dependent  plasticity  (creep),  the  instability  could  be  sup- 
pressed by  using  a pattern  of  varying  time  increments. 
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1.  Introduction 

There  is  a broad  need  for  understanding  the  fracture  of  structural 
alloys  at  stress  levels  approaching  general  yield.  Many  structures,  pa- 
ticularly  in  military  and  transportation  fields,  are  subjected  to  such 
loads  in  service.  Under  repeated  loads,  giving  low  cycle  fatigue,  one 
may  need  to  predict  the  direction  of  crack  growth  as  well  as  the  rate. 

In  some  cases,  fracture  is  intended,  as  in  metal-working  operations  or  in 
fragmentation.  In  other  cases,  one  would  like  to  predict  failure  in 
large,  elastic  structures  from  small,  nearly  plastic  test  specimens. 
Finally,  one  would  like  to  be  able  to  tell  metallurgists  how  changes  in 
alloys  might  improve  ductility  against  crack  initiation,  or  toughness 
against  crack  propagation.  This  report  summarizes  some  work  directed 
toward  such  ends. 

We  first  consider  the  different  mechanisms  of  fracture  in  typical 
structural  alloys  under  various  notch  configurations.  Ideally,  from 
such  information,  fracture  behavior  could  be  derived  by  mechanics.  Cur- 
rently, service  behavior  must  be  predicted  from  test  specimens,  so  we 
look  at  the  range  of  validity  of  the  J Integral.  We  correlate  differ- 
ent fracture  mechanisms  with  the  different  stress  states  that  are  found 
when  the  criterion  is  not  met. 

For  crack  growth  by  hole  coalescence  we  present  an  approximate 
analysis.  For  the  cleavage  fracture  mode  we  show  by  s^ome  experiments 
on  cracks  and  notches  at  defective  welds  in  steel  that  present  test  con- 
figurations may  not  be  conservative,  indicating  the  need  for  a wider 
range  of  test  specimen  configurations. 


Finally,  we  report  on  the  prospects  for  development  of  computer 
programs  useful  for  further  numerical  work  either  for  structures  to  de- 
termine the  stress  and  strain  in  the  neighborhood  of  the  fracture  pro- 
cess zone,  or  for  microstructures  to  estimate  crack  formation  and  growth. 
A condition  for  stable  explicit  integration  of  creep  equations  is  sug- 
gested. 

2.  Fracture  Mechanisms 

Fracture  initiation  in  metals  under  monotonic  load  has  been  in- 
vestigated over  a number  of  years.  The  fracture  site  in  the  center  of 
a necked  specimen  was  observed  by  Ludwlk  (1926)  and  an  approximate 
stress  analysis  carried  out  by  Bridgman  (1952).  Cleavage  fracture  in 
steel  at  low  temperatures  was  recognized  early.  The  mechanism  of  hole 
growth  from  inclusions  was  identified  by  Tipper  (1949)  although  it  was 
the  electron  micrographs  of  Plateau,  Henry,  and  Crussard  (1956)  that 
brought  a wide  understanding.  The  mlcromechanlcs  of  growth  of  non- 
interacting  holes  was  developed  by  McClintock  (1968) , Rice  and  Tracey 
(1969), and  Tracey  (1971),  but  tends  to  over-estimate  the  ductility.  Pre- 
liminary estimates  of  the  limitation  of  ductility  by  localized  flow 
between  holes  were  made  by  Nagpal,  et  al.  (1973),  but  better  analysis 
is  needed  for  the  beginning  of  concentrated  flow  by  hole  Interaction  in 
metals  with  a history  of  slight  strain-hardening. 

For  crack  growth,  several  mechanisms  have  been  reported. 

a)  A sharp  crack  growing  either  by  continuous  blunting  (Rice  and 
Johnson,  1970),  (Joyce,  1968),  by  alternating  sliding  off  as  a sharp  V 
(Joyce,  1968),  (Elliott  and  Stuart,  1968),  (Hayden  and  Floreen,  (1969), 
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or  by  several  such  wedges. 

b)  A crack  growing  by  successive  Interactions  with  a single  hole, 
as  suggested  by  Clpolla  {1973) . 

c)  Individual  voids  so  numerous  as  to  constitute  a continuum. 

d)  The  localization  of  flow  In  a dilating  continuum  (Berg,  1970), 
followed  by  the  joining  up  of  the  Individual  voids  within  the  band 
(Nagpal,  et  al,,  1973)  giving.  In  effect,  traction-displacement  boundary 
conditions  on  the  surrounding  material. 

Means  of  handling  each  of  these  crack  formation  and  growth  mecha- 
nisms must  be  Included  In  any  general  numerical  program  by  allowing  for 
concentrated  flow  elements,  arrays  of  holes,  dllatlonal  plasticity,  gen- 
eral compliant  and  Interactive  boundary  conditions.  Experiments  are 
still  required,  however,  to  decide  which  mechanisms  are  most  Important 
and  therefore  which  should  be  analyzed  first. 

3.  Valid  J-Integral  Specimens  for  Crack  Instability 

(Summary  of  Joyce  and  McClintock  (1976),  with  added  notes  on 
fractography) , 

Tests  on  a variety  of  alloys,  with  the  specimen  configuration  of 
Fig.  3,1,  led  to  a criterion  for  valid  J-correlation  of  crack  instability. 

Assume  an  Ilyushin  stress  strain  curve  with  exponent  n , stress  at  unit 
strain  , fitted  to  the  total  stress-strain  curve  at  the  yield  point 
and  the  fracture  point: 

a “ (e)"  . (3.1) 

I 

The  correlation  Is  based  on  the  Idea  that  a certain  stress  ratio  Is  { 


j 
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required  to  develop  the  Hutchinson-Rici'-Rosengren  (HRR)  non-l incar  strcs 
field,  within  a region  small  compared  to  the  specimen  ligament  dimension 
b but  large  compared  to  the  crack  tip  opening  displacement,  discussed 
further  In  Section  5, 

CTOD  = 2 . (3.2) 


The  required  radius 

HRR 

in  terms  of  a stress  ratio 

^ ’^HRR 


for  a valid 
c by 

_ "Ic  ( 


HRR 


singularity  is  then  given 

(n+1) /n 

(3.3) 


With  a value  of  c * 1.1,  the  criterion  successfully  distinguished 
between  specimens  with  valid  and  Invalid  J-correlations,  as  judged  by 
the  same  and  the  same  fractographic  features  for  both  the  single 

and  the  double  grooved  specimens.  Similar  success  was  found  for  a range 
of  steels,  except  that  for  the  higher  strength  steel,  the  factor  c must 
be  taken  to  be  1.05.  While  the  difference  appears  small,  the  low  expo- 
nent n makes  this  difference  in  constants  c significant.  The  seml- 
emplrical  nature  of  this  correlation,  and  the  fact  that  it  does  not  take 
into  account  differing  states  of  stress  and  strain  in  the  corresponding 
non-hardening  solutions  , mean  that  further  experiments  and  also  numerical 

calculations  of  the  stress  and  strain  fields  for  strain  hardening  expo- 

\ 

nents  of  the  order  of  0.1  would  be  very  desirable. 

A.  Fractographic  Observations 

Specimens  of  202A-T4  and  707S-T651  aluminum  alloys  with  ligaments 
as  small  as  1.91  mm  gave  apparently  valid  values  of  0.026  and 
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0.036  MJ/m  . Values  of  K,  calculated  from  these  values  were 

Ic  Ic 

found  to  be  larger  by  30  to  50%  than  handbook  values  for  these  ma- 

terials. This  difference  might  have  resulted  from  the  different  forming 
processes  of  the  specimens.  Any  difference  due  to  crack  initiation  oc- 
curring before  maximum  load  in  these  materials  is  likely  to  be  less  than 
20%,  judging  from  fatigued-loaded-fatigued  specimens.  The  initial  frac- 
ture topography  of  the  single  and  double  grooved  specimens  was  identical. 
The  crack  tip  blunting  and  the  average  dimple  diameter  are  on  the  order 
of  the  CTOD,  in  agreement  with  Rice  (1973).  Fracture  was  by  void  growth 
and  coalescence  and  occurred  before  the  limit  load  was  reached. 

The  fatigued- loaded-fatigued  specimens  of  2024-T4,  7075-T651,  and 

6061-T6  showed  that  at  most  one  hole  was  opening  in  front  of  the  crack 

\ 

tip.  No  bands  of  localized  flow  were  observed  which  had  grown  in  advance 
of  the  crack  tip. 

In  6061-T6  aluminum  and  cold-finished  1117  steel,  for  which  valid 
results  were  not  obtained,  the  initial  fracture  surface  topographies 
for  single  and  double  grooved  specimens  were  noticeably  different.  The 
double  grooved  specimens  showed  delamination,  a sharp  step  at  the  fatigue- 
fracture  transition,  and  then  an  abrupt  transition  to  hole  growth  fracture. 
The  step  height  was  on  the  order  of  the  CTOD,  which  in  turn  is  equal  to 
approximately  five  dimple  diameters.  On  the  single  grooved  specimens,  no 
step  was  observed  at  the  fatigue-fracture  transition.  Instead,  a combi- 
nation of  hole  growth  and  sliding  off  is  observed.  A transition  to  pure 
hole  growth  then  occurs  about  20  dimple  diameters  from  the  root  of  the 
fatigue  crack. 

In  a hardened  hlgh-purlty  aluminum  alloy,  valid  J tests  were  not 
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found.  Initial  fracture  topographies  wore  very  similar  for  single  and 
double  grooved  specimens.  In  specimens  with  cracks  normal  to  the  width 
direction  W , growing  in  the  thickness  direction  T (WT  orientation) , and 
also  in  ones  with  the  RT  orientation,  fracture  occurred  by  delamlnatlon 
consisting  of  hole  growth  in  the  plane  normal  to  the  direction.  The 
average  dimple  spacing  was  0.8  vim  and  inclusions  of  approximately 
0.15  ym  diameter  were  visible  in  the  center  of  many  of  the  dimples.  The 
TW  specimens  do  not  show  a hole  growth  fracture  mechanism,  which  is  a 
surprise  since  the  fracture  occurred  in  the  same  plane  which  delaminated 
by  hole  growth  in  the  WT  and  RT  specimens.  Instead,  the  fracture  sur- 
face consisted  of  steps  with  a wavy,  featureless  appearance. 


5.  Predicting  Crack  Toughness  from  Microstructure 
(Excerpts  from  McClintock,  1977) 


Take  as  a measure  of  the  crack  toughness  of  a material  the  critical 
plastic  zone  size  for  instability,  r^^  given  in  terms  of  the  critical 
stress  intensity  factor  and  the  yield  strength  Y for  a nonharden- 

ing material  by  Levy,  et  al  (1971) : 


r 


Yc 


(5.1) 


Most  alloys  whose  fracture  is  of  concern  are  so  lightly  strainhardening 

(n  z 0.1)  that  Eq.  5.1  is  a good  approximation. 

It  is  the  critical  value  of  this  plastic  zone  size.  r..  . not  K_ 

— Yc  Ic- 

Itself,  that  should  be  tabulated  as  a material  property  characterizing 

crack  toughness  because  its  dimensions  are  those  of  length  and  are  much 

-3/2 

more  easy  to  visualize  than  FL  . In  particular : 
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1.  The  designer  can  directly  compare  r^^  with  the  dimensions  of 
his  part.  If  the  calculated  r^^  is  greater  than  the  ligament  dimension, 
the  part  will  carry  its  fully  plastic  limit  load  before  fracture  (although 
fracture  may  still  reduce  the  extension  of  the  part). 

2.  r^^  gives  a rough  measure  of  the  maximum  allowable  half-length 
of  a crack  that  can  withstand  working  stresses  up  to  a large  fraction  of 
the  yield  strength.  Specifically,  up  to  about  half  the  yield  strength, 
the  relation  for  Kj.  in  terms  of  crack  half  length  c and  applied  stress 
at  "Infinity"  ^oo  t 


, 


(5.2) 


and  Eq,  5.1  combine  to  give 

c 

c 

For  instance,  for  a maraglng  steel  with  a reasonably  tough-sounding  value 
of  = 90  (83,000  Ib/in.^^^)  and  a yield  strength  of  1806  MN/m^ 

(262,000  psl),  the  critical  plastic-zone  size  is  only  0.39  mm.  Cracks 
down  to  about  this  size  would  have  to  be  found  and  eliminated,  unless  one 
adopts  fail-safe  design  or  is  willing  to  pay  off  on  failures. 

3.  The  plane  strain  r^^  can  be  compared  directly  to  the  specimen 
thickness.  The  plane-strain  assumption  is  valid  if  the  part  is  thicker 
than  about  6 times  the  plane-strain  critical  plastic-zone  radius. 

4.  If  r^^  approaches  the  ligament  dimension,  nonlinear  elastic 
fracture  mechanics  (the  J-integral)  must  be  used. 


5.  For  nonhardening  materials,  the  f lank-to-flank  crack-tip  opening 
displacement  can  be  readily  pictures  In  terms  of  the  plastic-zone  radius  by 


CTOD 


2.7  Y/E  . 


(5.4) 


_<j_ 


5.1  Displacement  criterion  for  crack  initiation 

Now  return  to  estimating  the  fracture  toughness  at  initial  crack 
growth,  r^j^  , from  the  inclusion  half  spacing  £ assuming 

CTOD^  = 2£  . (5.5) 

For  a nonhardening  material,  Eq.  5.4  gives 


r.^^  = 0.37  (2i)  (Y/E)  . (5.6) 

For  a hardening  material,  the  near  tip  stress,  strain,  and  dis- 
placement fields  can  be  given  in  terms  of  the  J-integral,  which  can  be 
thought  of  as  a scalar,  like  . These  fields  hold  for  a nonlinear, 
elastic,  incompressible  stress-strain  relation  of  the  form 


o 


(5.7) 


If  Eq.  5.7  represents  the  total  strain,  the  nonlinear-elastic  or  de- 
formation theory  plasticity  gives  the  same  result  as  the  more  physical 
Incremental  plasticity.  It  is  not  valid  for  drack  growth,  wherein  the 
ratios  of  strain  increments  change  as  the  crack  passes  by.  For  most 
alloys,  the  exponent  n is  of  the  order  of  0.05  to  0.25.  The  con- 
stant can  be  thought  of  as  the  equivalent  stress  at  unit  strain. 

The  nonlinear  stress  and  strain  fields  also  depend  on  a numerical  co- 
efficient  I^  that  varies  slowly  with  n within  2 percent  of  (Joyce  and 
McClintock,  1976) 


I 

n 


10.3  /0.13  + n - 4.8n  . 


(5.8) 
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In  terms  of  these  quantities  and  normalized  stress,  strain,  and  displace- 
ment functions  which  depend  only  on  angle,  the  Hutchinson  (1968a, b)  Rice 
and  Rosengren  (1968)  (HRR)  fields  are 


n/(n+l)^ 


aij(e) 


^j(r,e)  “ 


,l/(n+l) 


0,1  r 
1 n 


ey(e) 


Uj^(r,0) 


1/  (n+1) 


Uj^(0) 


(5.9) 


At  the  flank  of  a crack,  the  normalized  displacement  function  u^(180°) 
Is  unity.  The  desired  flank- to- flank  opening,  where  r = u , Is  then 

CTOD  - • (5.10) 

n 1 

Equation  5.10  gives  perhaps  the  most  direct  physical  Insight  Into  the 
meaning  of  J.  Since  the  deformation  theory  stress  strain  relation  of  Eq. 
5.6  Is  not  at  all  valid  for  crack  growth,  the  Interpretation  of  J as  a 
strain-energy  release  rate  Is  not  physically  very  useful,  even  though  It 
Is  provided  the  Insight  for  the  original  derivation. 

When  the  plastic  zone  Is  extending  nearly  all  across  the  specimen, 

J must  be  evaluated  numerically  for  the  particular  structure  and  loading. 
When  there  Is  a surrounding  linear  elastic  stress  singularity,  J can 


be  found  from  K : 
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J = Kj  (1  - v^)/E  . (5.11) 

In  terms  of  the  nonhardening  plastic-zone  radius, 

J = 2Tr  r^  (1  - v^)/E  . (5.12) 

Finally,  we  can  relate  the  critical  plastic-zone  size  for  crack  initiation 
to  the  spacing  between  hole  nuclei  by  combining  Eqs.  5.5,  5.10,  and 
5.12: 

I Ea 

^ r . with  V = 0.5  . (5.13) 

2tt  - v2) 

For  a nonhardening  material  this  reduces  to  r^j^  = 0.39  (2Jl)E/Y,  compared 
to  0.37(2il)E/Y  obtained  from  Eq.  5.6.  For  a hardening  material  with 
n = 0.1,  Y/E  = 0.004,  'i/a^  = 0.174»  r^^  = 0.822(2il)E/Y,  indicating  a 

rather  strong  effect  of  strain  hardening. 

5.3  Crack  Instability 

There  is  not  much  prospect  of  a closed  form  analytical  solution 
for  a growing  Mode  I crack,  taking  into  account  geometry  changes,  the  in- 
cremental flow  rule,  strain  hardening,  residual  stresses,  and  the  strained 
material  left  behind  the  growing  crack.  Perhaps  some  Insight  can  be  ob- 
tained from  the  Mode  III  solutions  (shear  parallel  to  the  leading  edge  of 
the  crack).  The  most  unrealistic  feature  is  that  in  Mode  III  the  plastic 
zone  extends  mostly  in  front  of  the  crack,  so  the  crack  advances  into  pre- 
strained material,  whereas  for  initial  Mode  I loading  the  plastic  zone 
runs  off  to  the  side.  This  discrepancy  is  likely  to  be  reduced  by  strain 
hardening,  residual  stress,  and  possible  zigzagging  of  a Mode  I crack. 


With  these  caveats,  the  Mode  I analogue  of  Mode  111  instability  is 
given  in  terns  of  the  fracture  strain  (perhaps  including  any  needed 

to  initiate  holes  and  only  the  further  strain  to  the  point  of  localiza- 
tion), McClintock  (1971): 


r 


Yc 


21  exp 


/(Y/E)  + 1 - 


(5.14) 


5,4  Strain  criterion  for  initiation.  The  analogue  of  Eq.  5.14  for 
initiation  is: 


r^j  = 21  (e^P/(Y/E)  + 1)  . (5.15) 

5,5  Comparison  with  data.  The  estimates  for  Initiation  by  dis- 
placement (Eq.  5.13)  and  strain  (Eq.  5.15)  and  for  instability  (Eq.  5.14) 
are  given  in  Table  2.5  for  the  three  aluminum  alloys  and  the  1117  cold- 
finished  steel  studied  by  Joyce  (1974) . The  hole  growth  ratio  was 
taken  roughly  from  his  scanning  electron  micrographs.  The  fracture  strain 
was  assumed  to  consist  entirely  of  that  required  to  grow  the  holes  to  co- 
alescence. An  approximate  equation,  exact  for  zero  or  linear  hardening 
and  equl-axial  transverse  strain  and  axial  shortening  of  cylindrical  holes 

has  been  obtained  by  McClintock  (1968) . The  rate  of  change  of  semi-minor 

00 

axis  a of  holes  spaced  2x.^  apart,  subject  to  applied  stress  of 
00 

and  a|^  and  internal  pressure  p , in  a material  with  an  equivalent  flow 

.00 

strength  Y undergoing  an  equivalent  strain  increment  of  e is 

d In(a/£  )/de"“  - slnhf  (1-n)  (a” -t-a*  + 2p)/(2Y/»'3)  ]/(l-n)  . (5.16) 

a o D 


Note  that  in  the  limit  as  the  strain-hardening  exponent  n approaches  0 


the  holt'  growth  rate  ht'eomes  exponential  wllli  the  applied  Htrean,  where- 
as for  linear  hardening,  n = 1,  the  hole  growth  rate  per  unit  applied- 
strain  increment  varies  linearly  with  the  applied  stress.  Equation  5.16, 
for  the  stress  ratio  (Og”  + )(2Y//5)  of  3/2  + tt  , expected  in  front 

of  a crack,  was  used  to  obtain  the  fracture  strain  needed  In  Eqs.  5.14 
and  5.15  at  the  observed  hole  growth  ratio. 


Table  5.1  Initial  and  Critical  Plastic  Zone  Radii 


!Alloy 

Y, 

MN/m^ 

E, 

MN/m^ 

n 

2£, 

TTTm 

Initiation,  f^j^ 

Instability,  f^^  * 

(Jl  /a) 

Strain  Dlsplac. 

Eq.5.13  Eq.5.15 
nun  non 

Strain 

Eq.5.14 

mm 

Observed 

nun 

i7065-T6 

500 

72,000 

0.10 

0.012 

3 

1.42 

0.064 

0.1 

1.34 

] 

1 

0.016 

4 

1.89 

0.104 

0.189 

1 

!202A-T4 

390 

72,000 

0.12 

0.012 

3 

1.82 

0.084 

0.163 

3.39 

! 

0.040 

3 

6.07 

0.281 

0.544 

6061-T6 

281 

72,000 

0.09 

0.012 

6 

2.53 

0.159 

0.689 

11. 

1117 

510 

200,000 

0.07- 

0.004 

5 

1.89- 

0.332- 

2.14 

>25. 

steel 

0.11 

0.020 

6.44 

0.377 

3.13 

t 

For  the  stronger  alloys  of  Table  1,  the  displacement-based  plastic 
zone  size  for  Initiation,  r^^  , is  comparable  to  the  critical  value,  as 
has  been  noted  by  others.  The  displacement-based  criterion  fails,  however, 
to  predict  the  very  large  Increases  in  plastic  zone  radii  in  alloys  with 
low  yield  strain  Y/E  . 

The  strain-based  criteria  show  relatively  little  difference  between 
initiation  and  instability  for  the  strongest  alloys,  but  much  more  for  the 
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weaker  ones,  in  line  with  experience.  The  theoretical  level  of  ty 
Is  too  low  for  all  the  strain  criteria.  Possible  reasons  include  crack 
blunting  with  loss  of  trlaxiallty,  lack  of  strain  hardening  In  Eq.  5.14, 
and  advance  of  the  crack  Into  unstrained  material.  Detailed  numerical 
work  appears  necessary,  in  addition  to  any  possible  analytical  results. 

6.  A Grooved-Strip  Test  for  Weld  Ductility 

6.1  Introduction 

The  question  examined  here  Is  whether  an  Improved  test  specimen 
can  provide  a better  measure  than  the  Charpy  V-notch  for  the  practical 
performance  of  plates  and  welded  structures.  Because  It  would  be  de- 
sirable for  welded  structures  to  become  fully  plastic  before  fracture, 
even  in  the  presence  of  a crack,  and  because  analytical  and  numerical 
methods  for  determining  the  ductility  of  a large  plate  in  terms  of  the 
Charpy  behavior,  or  critical  stress  Intensity  factor,  for  example,  are 
not  known  for  the  fully  plastic  case,  Marcollni  (1975)  carried  out  a 
test  on  specimens  containing  defects  as  near  to  full  scale  as  practical. 

6.2  Material  and  welding. 

The  chemical  composition  of  the  A-36  base  metal  was  0.21%  C, 

0.69%  Mn,  0.04%  Si,  0.015%  P,  and  0.018%  S.  Properties  In  13  mm  plate 
were  279  MN/m^  (40,400  psi)  yield,  439  MN/m^  (63,600  psi)  tensile 
strength,  and  27%  elongation  In  203  mm  (8  In.). 

Two  manual  submerged  arc-weld  passes  were  made  with  a 2 mm  (5/64  in. 
diameter  L-60  electrode/860  flux  combination  (Llncolnweld  R,  Classifica- 
tion F62-EL12).  Tests  required  by  specifications  AWS  A5. 17-69.  ASME 
SFA  5.17  showed  the  chemistry  and  mechanical  properties  of  the  deposited 
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wold  motal,  aocordlri);  iv»  tho  tost  co«  i i fli-aio  lor  this  ( i rmU’ 
combination,  to  be  0.0772  C,  0.882  Mn,  0.222  Si,  0.0192  P , and  0.0152  S , 
with  a yield  of  405  MN/m^  (58,800  psi),  a tensile  strength  of  484  MN/m^ 
(70,200  psi),  an  elongation  in  51  mm  of  292,  and  Charpy  V-notch  impact 
properties  of  85  Nm  (63  ft  lb)  at  244  K (-20  F) . 

For  the  largest  practical  test  with  available  equipment,  strips  of 
plate  12  mm  thick  by  76  mm  wide  were  welded.  Surface  cracks  were  simu- 
lated by  machining  grooves  2 mm  deep  for  9.5  mm  with  a 102  mm  dla.  cutter, 
in  the  face  of  the  strip  and  parallel  to  the  weld,  at  various  locations 
base  metal,  weld  metal  and  heat-affected  zone  (HAZ) . The  strips  were 
then  subjected  to  alternating  bending  in  an  SFl  fatigue  testing  machine 
of  ±340  Nm  (3000  in. -lb)  capacity.  Strips  of  this  size  could  be  pulled 
in  a 0.9  MN  (200,000  lb)  hydraulic  testing  machine.  Adequate  ductility 
is  indicated  by  a significant  reduction  of  area  in  the  minimum  section 
and  a hole  growth  mode  of  fracture  as  the  crack  progresses  out  of  the 
pre-cracked  seml-elllptlcal  zone  into  the  unnotched  edges  of  the  strip. 

6.3  Results 

The  principal  results  are  shown  in  Table  3.1.  Charpy  tests  showed 
a severe  degradation  of  behavior  in  the  heat-affected  zone,  where  the 
transition  temperature  at  least  was  294  K,  as  contrasted  to  255  to  263  K 
in  the  base  and  weld  metals.  For  the  singly  grooved  strip,  the  fracture 
appearance  transition  temperature  for  the  base  metal  and  the  heat  affected 
zone  were  much  more  nearly  comparable  (250-260  K) . The  transition  tem- 
peratures as  indicated  by  thinning  of  either  the  center  section  or  the 
unnotched  edge  section  were  very  nearly  the  same  as  the  fracture  appear- 
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ance  transition  temperature,  This  is  likely  because  all  three  tran- 
sitions were  concerned  with  a growing  crack.  With  doubly  grooved  spe- 
cimens, however,  (which  did  not  contain  a fatigue  crack  extension)  the 
fracture  appearance  transition  temperature  after  the  crack  began  to  run 
was  greater  than  299  K.  The  crack  was  ductile  in  traversing  the  thick- 
ness, but  changed  to  cleavage  as  it  spread  laterally.  Note  that  a speci- 
men as  small  as  a Charpy  specimen  would  have  indicated  complete  ductility 
under  these  circumstances.  It  would  erroneously  have  indicated  safe  be- 
havior. In  the  wider  strip  the  crack  apparently  became  unstable,  picked 
up  speed  and  changed  to  a cleavage  mode,  which  would  be  catastrophic  in 
a large  structure.  Similar  doubly-grooved  behavior  was  observed  for 
the  base  metal. 

6 .4  Conclusion 

Apparently,  the  Increased  trlaxiality  due  to  two  grooves  opposite 
each  other,  led  to  a 40  K rise  in  fracture  transition  temperatures  for 
A-36  steel.  Furthermore,  the  material  does  not  respond  to  these  grooves 
potential  defects)  in  the  same  way  that  the  Charpy  specimen  does.  The 
Charpy  specimen  may  be  either  over-  or  under-conservative,  depending  on 
the  conditions.  There  is,  therefore,  a need  for  Improved  tests  that  more 
» closely  simulate  service  conditions,  until  such  time  as  theory  and  numeri- 

cal methods  can  provide  quantitative  predictions.  Similar  results  would 
be  expected  in  fragmentation. 
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Table  6.1 

Fracture  data 

for  welded,  i^rooved 

strips  of  A-36  steel 

12.7  X 

76  mm, 

610  mm  long,  254  mm 

between  grips 

No.  Total 
crack 
depth 

Temp. 

Load 

rate 

A TS 
Sr 

Center  Edge  Crack  Fracture 

thinning  thinning  opening  mechanism 
angle 

mm 

K 

-1 

s 

% X 

0 

Base  metal 


Base 

-1 

8.7 

247 

2.91 

E-3 

1.6 

6.1 

0 

cleavage 

-5 

6.8 

248 

0.95 

ft 

3.9 

4.4 

4 

mixed 

-6 

5.9 

250 

1.39 

If 

8.4 

18. 

15 

hole  growth 

-4 

7.7 

256 

1.54 

ft 

8.7 

19. 

24 

hole  growth 

-3 

8.7 

256 

0.21 

If 

8.7 

17. 

22 

hole  growth 

-2 

7.1 

265 

0.92 

If 

6.8 

5.6 

5.4 

mixed 

-7 

4.1 

299 

1.50 

II 

1.2 

4.4 

mixed 

(Double  groove,  2 mm  each  side,  no  fatigue  cracking.) 


Heat-affected  zone 


HAZ 


-4 

8.0 

247 

2.62 

E-3 

0.0 

3.5 

0 

cleavage 

-7 

6.4 

249 

0.82 

II 

0.0 

1.7 

0 

cleavage 

-6 

6.2 

258 

0.70 

II 

9.4 

14.1 

18 

hole  growth  ^ 

-5 

6.9 

263 

0.98 

If 

8.3 

10.2 

25 

-2 

5.8 

275- 

286  0.1-0 

.16" 

8.4 

12.9 

19 

-3 

3.2 

280- 

•286  0.1-0 

.16" 

6.2 

6.6 

10 

-1 

5.9 

286- 

292  0.1-0 

.16" 

9.4 

13.1 

19 

-9 

5.5 

299 

18.9  " 

7.2 

10.1 

10 

1 

-8 

3.7 

303 

12.2  " 

9.8 

18.2 

17 

i 

-10 

4.1 

299 

1.3 

0.0 

3.5 

- 

mixed 

(Double  groove,  as 

above) 

4 
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7.  Numerical  Studies 

Joyce  (1974)  developed  a numerical  method  for  plastic  plane  strain 
analysis  based  on  boundary  Integral,  stress  relaxation  techniques.  The 
method  was  extended  to  creep  by  Eblsu.  The  method  Is  based  on  an  elastic 
boundary  Integral  program  Implementing  the  following  scheme: 

a)  Boundary  conditions  can  be  of  the  most  general  linear  type, 

It 

giving  a number  of  loading  parameters  P as  linear  functions  of  the  com- 
ponents of  traction  (force  per  unit  surface  area,  t^)  and  displacement  u^ 
of  each  boundary  segment  y . 

pk  ^ p ky  ^ p ky  . (7.1) 

t u 

(For  example,  traction  boundary  conditions  are  met  by  setting  P^.  = 1 , 

ali  other  terms  zero) . These  general  boundary  conditions  allow  treating 
the  fracture  mechanisms  of  Section  2.  It  would  be  very  desirable  to  ex- 
tend the  capabilities  of  finite  element  programs  to  Include  them.  (For 
the  boundary  Integral  method.  In  order  to  satisfy  moment  equilibrium.  It 
Is  necessary  to  Include  moments,  normally  zero,  applied  to  each  segment 
which  are.  In  effect,  gradients  of  the  normal  traction,  t^^  . (Shear 
gradients  of  traction  and  gradients  of  displacement  are  Included  for  com- 
pleteness and  Improved  accuracy,  although  for  simplicity  they  will  not  be 
Included  In  the  following  discussion). 

b)  With  the  body  modeled  as  being  embedded  In  an  Infinite  matrix, 
relative  displacements  are  Inserted  between  the  body  and  the  matrix.  At 
either  end  of  any  given  segment  m there  Is  then  a pair  of  edge  disloca- 
tions, called  darts,  D™  . The  known  stress  fields  of  dislocations  give 
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the  tractions  t“  and  displacements  on  the  at' tooted  soktnionts  n 

due  to  the  dislocated  sei5inonts  ro  . Slimming  over  all  si'p.nu'nts. 

jP  = tP™  d”  , D™  . (7.2) 

c)  Substitution  of  Eqs.  7.2  into  7.1  gives  an  equation  for  the 
unknown  darts  in  terms  of  the  parameters  of  the  boundary  conditions, 
which  can  then  be  solved  for  the  darts: 

pk  , rpkp^pm  ^ p ky  ^yn^  ^m  _ 

d)  With  the  darts  known,  the  stresses  at  any  interior  points  can 
be  found. 

e)  Plasticity  is  modeled  by  subdividing  the  interior  into  elements. 
Within  the  elements,  strain  is  modeled  by  dislocations  migrating  to  the 
element  walls.  The  resulting  dislocation  arrays  reduce  the  stresses  and 
modify  the  boundary  tractions  and  displacements.  Their  effects  are  added 
to  Eqs.  7.2  and  the  process  is  repeated. 

The  program  was  developed  to  a working  state,  but  it  turned  out 

that  Instabilities  developed  after  long  times.  These  instabilities  also 

4 

appear  in  explicit  integration  of  rate-dependent  plasticity  (creep)  with 
the  finite  element  method.  Eblsu,  in  unpublished  work  supported  by  the 
project  showed  that  short  time  increments  re-stablllzed  the  calculations, 
allowing  a few  long  time  Intervals  before  instability  reappeared.  In  rate- 
independent  plasticity,  a subsequent  one-dimensional  analysis  indicates 
that  the  method  is  Inherently  unstable,  although  the  instability  shows  up 
slowly  enough  so  that  some  problems  can  be  solved. 


While  the  boundary  integral  method  developed  here  requires  solving 
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fewer  simultaneous  equations  than  the  finite  element  method,  the  matrices 
are  fully  dense,  not  sparse,  and  very  large  numbers  of  Influence  coef- 
ficients are  required.  It  appears  to  be  useful  primarily  for  small  plas- 
tic zones  around  cracks  In  large  bodies.  For  Instance,  It  has  been  suc- 
cessfully applied  to  the  cracks  In  railroad  rails,  under  complex  histories 
of  compression  and  shear  (McCllntock,  1977). 

8.  Conclusions 

1.  Singly  and  doubly  grooved,  plane  strain  tensile  specimens  cut 
from  25  mm  D bar  stock  were  used  to  test  the  validity  of  the  Hutchinson- 
Rice-Rosengren  (HRR)  J- integral  correlation.  A criterion  based  on  a de- 
sired stress  ratio  in  the  HRR  field  Indicates  valid  results  if  the  stress 
at  the  half-ligament  dimension  Is  no  more  than  1.05  to  1.1  times  the 
yield  strength.  For  smaller  specimens,  the  J integral  varied  linearly 
with  ligament  size  and  was  surprisingly  similar  for  single  and  double 
grooved  specimens. 

2.  In  the  valid  J region  for  high  strength  aluminum  alloys,  crack 
Initiation  occurred  at  a crack  tip  opening  of  the  order  of  one  hole 
spacing.  Crack  growth  occurred  by  successive  interaction  of  one  hole  at 

a time  with  the  crack.  In  specimens  of  6061-T6  aluminum  and  cold  finished 
1117  steel, with  Invalid  J criteria,  single  and  double  grooved  specimens 
failed  by  different  mechanisms.  The  double  grooved  specimen  showed  de- 
lamlnatlon,  a sharp  step  at  the  point  of  fracture  Initiation,  and  then  an 
abrupt  transition  to  hole  growth.  For  6061  T6,  the  step  height  was  of  the 
order  of  five  hole  diameters.  Only  on  the  single  grooved  specimens  of 
these  alloys  was  sliding  off  observed.  The  transition  to  hole  growth  oc- 
curred after  20  hole  diameters  of  CTOD. 
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3.  Plates  of  a high  purity  aluminum  zinc  alloy  exhibited  very  an- 
isotropic fracture  behavior,  sometimes  with  a wavy  featureless  appearance 
rather  than  hole  growth. 

4.  A theoretical  prediction  of  the  critical  plastic  zone  size  from 
fractographic  observations  shows  that  the  displacement  concept  failed  to 
indicate  the  very  large  Increases  in  plastic  zone  in  alloys  with  low  Y/E. 
A criterion  based  on  a critical  strain  and  an  analogy  with  Mode  III  in- 
dicates the  correct  effect  of  Y/E  but  falls  below  the  observed  critical 
radius. 

5.  In  regard  to  the  transition  temperature  for  welded  low  carbon 
steel,  the  effect  of  plastic  constraint  Induced  by  double  grooved  speci- 
mens can  cause  a 40®C  rise  in  fracture  transition  temperature  above  that 
observed  by  the  Charpy  specimen.  Thus  under  some  conditions  the  Charpy 
test  may  be  grossly  unconservative. 

6.  Studies  of  numerical  methods  based  on  a boundary  Integral 
method  and  dislocation  modeling  of  plasticity  within  elements  showed 

an  instability  which  could  be  suppressed  in  creep,  but  required  enough 
computation  with  influence  coefficients  so  that  it  does  not  appear  com- 
petitive with  finite  element  methods,  except  for  small  plastic  zones 
around  cracks  in  large  bodies. 
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